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SUMMARY 

The performance of a class of remotely piloted, microwave-powered, high- 
altitude airplane platforms was studied. Each cycle of the flight profile con- 
sists of climb while the vehicle is tracked and powered by a microwave beam, 
followed by gliding flight back to a minimum altitude. Parameter variations 
were used to define the effects of changes in the characteristics of the air- 
plane aerodynamics, the power-transmission systems, the propulsion system, and 
winds. 

Results show that wind effects limit the reduction of wing loading and 
increase of lift coefficient, two effective ways to obtain longer range and 
endurance for each flight cycle. Calculated climb performance showed strong 
sensitivity to some power and propulsion parameters. A simplified method of 
computing gliding endurance was developed. 


INTRODUCTION 

Remotely piloted vehicles operating at high altitude have been proposed to 
perform communication or observation tasks for various regions of the Earth's 
surface (refs. 1 and 2). A remote power supply, such as solar radiation or a 
microwave beam from a ground station, could give endurance limited only by 
systems reliability. Applications for such high-altitude aircraft platforms 
include mapping, resource monitoring, relaying communications, and conducting 
other tasks currently performed by satellites or manned aircraft. 

Long-endurance aerial platforms offer advantages over alternate systems. 
Endurance of a manned aircraft is limited by fuel storage and human fatigue. 
Furthermore, the payload and equipment must include provisions for the pilot 
and his environmental control system. These factors adversely affect cost and 
complexity. A geosynchronous satellite has long endurance; however, it also 
has high cost, poor resolution for observation tasks, and constraints for com- 
munications tasks because of the extreme range. A satellite operating In a low 
orbit passes only infrequently and briefly over a given region. Compared to a 
platform in the upper atmosphere, a low-orbit satellite requires observation 
systems to have resolutions at least six times as great for equivalent quality 
of results. i 

Several types of high-altitude aircraft platforms have been proposed. 
Lighter-than-air concepts have been considered (ref. 3). Some of the difficul- 
ties of operating these vehicles at altitudes of 18 km (60,000 ft) and above 
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[flf® environment. The airships would have to generate 

lift at air densities less than one-tenth that of sea level (ref. 41 and 
according to reference 5, have the capability to fly against winds of up to 

m/s (100 kt). Airplane configurations using solar power have 
been discussed in references 6, 7, and 8. The study of reference 8 concluder 

fn ? ih/r?2?\ technology and extremely low wing loadings (down to 15 Pa 

,)) would be required. Even so, a solar-powered configuration would 
be constrained to operate at low latitudes to obtain enough daylight hours and 
at locations and altitudes with modest wind velocities. 


Studies of the design 
airplane platforms (HAAP) 
HAAP configurations of all 
ating in the low-wind regi 
in the wing lower surface 
current power for electric 
were feasible based on the 
microwave technology (such 


and operation of microwave-powered high-altitude 
have been reported in references 3 and 9 to 11. The 
these reports were propeller-powered airplanes oper- 
on near 20 kilometers (66,000 ft) altitude. Rectennae 
receive microwave energy and convert it to direct- 
motors. These studies indicated that the designs 
assumption of some extrapolation of existing 
as that described in ref. 12). 


uAAD^^ predictions of cruise performance for the class of 

HAAP configurations that use a "linear" mode of flight. (These are described 
in the feasibility study of reference 10 and the design sensitivity study of 
reference 10.) In that mode, the flight profile consists of poSd climb near 

Station, followed by gliding flight that either returns the 
vehicle to the same microwave ground station or carries it to another station 

^tudy. Emphasis ?s placed on 

vehicle design and not power transmission. 

to Pa^ametric study define trends that apply 

to HAAP vehicles within a wide range of sizes and weights. Performance is 

of/HUude at the end of climb, exce™ eSerJJ 
Stored in batteries, range, and endurance for each cycle of the flight profile 

i Hcf variations in the aerodynamic ?hara^teJ- 

istics, power-transmission system, propulsion system, flight profile and 
winds. A minimum altitude of 18 kilometers (59,000 ft.) was Elected for all 
cases as a probable constraint due to wind. 

^ mi crowavc-powered airplane are sufficiently 
unconventional to require the development of a new computer' program for perfor- 
ce prediction. The program used in this study is documented in Appendix A. 


SYMBOLS 

Positive senses of some angles, axes and forces are presented in figure 1. 
A wing aspect ratio 

Ap propeller disk area, m2 

3 constant defined in equation (11) 
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wing span, m; also, constant defined in equation (11) 
drag coefficient, D/qs 
profile-drag coefficient 

lift coefficient, L/qs 
propeller-power coefficient, Pp/pn^Dp^ 

drag, N 

propeller diameter, m 
stored energy, J 
total energy received, J 

Oswald efficiency factor; also, base of natural logarithms 

acceleration of gravity, 9.80 m/sec2 at sea level 

altitude above sea level, km 

altitude at initiation of glide, km 

altitude at beam intercept point, km 

propeller advance ratio, V/nDp 

acceleration correction factor (see eq. (4)) 

microwave-beam intensity factor (see eq. (12)) 

wind-profile scale factor 

lift, N 

propeller rotational speed, revolutions/second 
maximum power available in beam, W 
power absorbed by the propellers, W 
pcwer available at rectenna, W 
dynamic pressure, 1/2 P V^, N/m^ 

reference value of radial distance from microwave ground station, km 

actual radial distance from microwave ground station, km 
2 

wing area, m^ 



T propeller thrust, N 

Tj degraded propeller thrust, N 

t elapsed time, s 

V true airspeed, m/s 

Vg equivalent airspeed, V yp/pg, m/s 

Vg ground speed, m/s 

Vtip propeller tip speed, m/s 

local horizontal wind speed, m/s 
W vehicle gross weight, N 

X horizontal range, km 

X3 horizontal distance between ground station and beam intercept point, km 

z dummy variable of integration, km (see eq. (B 3 )) 

Y flight path angle, deg 

^ increment of parameter 

n propeller efficiency factor 

6 microwave-beam/evaluation angle, deg 

M angle between wind vector and required ground-track vector 

P air density, kg/m^ 

pQ sea-level air density, 1.255 kg/m^ 

angle between airplane heading and required ground-track vector 
Subscripts 

c end of climb 

g end of glide 

max maximum 

t total for one climb and glide cycle 

A dot over a symbol denotes differentiation with respect to time. A bar 
over a symbol denotes an average value. 



CONCEPT DESCRIPTION 


The remotely piloted, microwave-powered HAAP of this study is based on the 
concept described in reference 9. Drawings of representative vehicles are 

shown in figure 2. A baseline configuration and microwave system are described 
in table I. 

The linear mode of flight, used in this study (and those of references 9, 
10, and 11) has a two-part cycle. The climb segment begins when a microwave 
beam begins to track the vehicle and transmit power to it. That power is used 
to climb and to store energy in batteries for use by the payload, guidance, and 
control systems. After power transmission terminates, the vehicle begins a 
long glide that either carries it to another ground station or back to the same 
station. 

The transmission of microwave energy is modeled largely with the assump- 
tions of reference 9. The multi-element retrodirective array or equivalent 
antenna (ref. Ip transmits a linearly polarized beam that is focused on the 
rectenna built into the HAAP. The two-dimensional tracking capability of the 
transmitter constrains the vehicle to fly in a predefined vertical plane over 
the ‘ground station. The sum of all range-related phenomenon is assumed to 
attenuate beam intensity as an inverse function of range. 

The conceptual design of the vehicle for this study is similar to a powered 
version of a high-performance sailplane. The payload fraction of 0.3 contains 
allocations for the observation or communications payload, batteries, and the 
guidance and control systems. The baseline configuration calls for high aero- 
dynamic efficiency to be achieved with high-aspect-ratio wings and extensive 
amounts of natural laminar flow. The wing-mounted rectenna uses linear polari- 
zation unless otherwise noted. Rectenna-received power below a minimum level 
for motor starting, or above the power capacity of the motor, is stored in 
batteries. Power in the range required for propulsion is used by high- 
efficiency electric motors to drive variable-pitch, constant-speed propellers. 
During both flight segments, the vehicle remains trimmed at one lift coeffi- 
cient. Although the propellers stop and fold streamwise when not in use to 
reduce drag, there is still a small increment in drag during gliding flight. 

A more detailed study of HAAP should consider criteria for stability, 
control, aeroelasticity, reliability, and other factors. The design illus- 
trated in figure 2(b) reflects some concern for reliability by minimizing the 
number of essential systems. Propulsion is provided by a single propeller. 
Aerodynamic control is achieved through control surfaces at the end of the tail 
booms; differential inputs of the horizontal surfaces produce wing twist. This 
design was examined briefly in an unpublished study which indicated that the 
two vehicles of figure 2 can have the same performance, control power, and 
weight. The generalized approach of this study is not constrained by choice of 
configuration. 
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ANALYSIS 


The evaluation of vehicle performance for a microwave-powered airplane 
requires mathematical modeling of vehicle motion, atmospheric effects (includ- 
ing wind), power transmission, and propulsion-system characteristics. The 
development of the equations used in the microwave-HAAP performance program of 
Appendix A is given in the following sections. 


Flight Mechanics 

Equations for force balance along the body axes can be developed with the 
conventions shown in figure 1(a). The associated assumptions are that: thrust 

and drag act along the same axis; airspeed is adjusted by changing flight-path 
angle to obtain the required lift; configuration lift coefficient remains 
constant; and excess thrust is used to climb. The resulting equations are: 

L - W cos Y = 0 ( 1 ) 

T-D-WsinY--|V = 0 (2) 


These equations can be modified to obtain the forms used in the performance 
program of Appendix A. First, thrust can be described in terms of propeller 
efficiency as: 

T ^ V'' (3) 

The term for the inertial acceleration force can be written as follows: 

W V = W dV ^ , W ^ V 3 ^.^ 
g g dh dt g dh 


For sufficiently small incremens of altitude, an acceleration correction factor 
Ka can be defined as: 


b ^ V dV ~ V AV 

^ g "dlT g aFI 


(4) 


Thus, equation (2) can be rewritten as 


n Pp/V - D - (1 + ka) W sin Y = 0 


(5) 


An expression for true airspeed, 


V = 



2 cos Y 
P Cl 


obtained from equation (1) may be written as: 

( 6 ) 


Except for the term cos Y, equivalent airspeed is simply the equilibrium 
airspeed for a given configuration at sea level: 


6 



( 7 ) 



The equations for the vehicle trajectory above a flat earth are based on 
the conventions shown in parts (b) and (c) of figure 1: 

V cos Ysin<J)- V„ sin y = 0 (8) 

X - V cosycos4> + cos y = 0 (9) 

h - V sinY = 0 (10) 

The use of these equations assumes that the airplane heading is automatically 
adjusted to compensate for the effects of wind. During climb, the resulting 
flight path must lie in the unique vertical plane swept out by the path of the 
microwave beam. 

Several parameters are functions of altitude. Air density is modeled on 
the geometric standard atmosphere of reference 4. The ratio of local density 
to sea-level density is calculated as 

p/Pq = (11) 

where the exponential coefficients hold constant over a typical altitude incre- 
ment of two kilometers. As in reference 10, it is assumed that L/D increases 
with altitude for the operating range of cruise altitudes because of the 
greater extent of laminar flow. The value of L/D is decremented (as a function 
of propeller size) for glide because of the drag of the folded propeller. 


Power Transmission 

The available power at the vehicle rectenna is assumed to be proportional 
to both range and the angular orientation of the rectenna surface with respect 
to the beam. Although the beam is considered to be focused, the effects of 
focusing precision and other factors are represented by a reciprocal 
relationship with range: 

Pj. oc (R/r)'^»' (12) 

where R is a reference radial range from the ground station, r is the true 
radial range, and kj. has a nominal value of 1.0. The power available is 
assumed to be proportional to the projected area of the rectenna that can be 
observed from the microwave ground station (eq. 122, vol. 2 of ref. 13). This 
can be approximated as 

P^ oc s sin (0+Y) ( 13 ) 

If both transmitter and rectenna use linear polarization, the power transfer 
can be conservatively approximated as a function of the phase angle between the 
two units (ref. 14 and eq. 25, vol. 1 of ref. 13): 
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(14) 


2 

Pp cc cos 4> 

If the rectenna or antenna have circular polarization, then the transmission 
efficiency drops by a factor of one half but remains unaffected by relative 
ground angle between the antenna and rectenna axes. 

All of these power-transmission relationships (expressions 12, 13, and 14) 
can be combined into one equation. It is convenient to describe power 
available for storage or propulsion in terms of power per unit weight as: 

-«'■ “ I (f )(7f 

where P/S is the maximum transmitted power per unit wing area available at the 
reference range R, W/S is wing loading, and both antenna and rectenna have 
linear polarization. As used in Appendix A, the equation, in addition, assumes 
an efficiency factor of 74 percent between the power reaching the rectenna 
surface and the power delivered to either the propeller shaft or the batteries. 


Propulsion 

Values of propulsive efficiency (eq. (3)) are determined from the tabulated 
values of reference 15 and are given as functions of advance ratio J and 
propeller-power coefficient Cp. These latter quantities can be determined as 
functions of both calculated and input parameters of the program of Appendix A: 



where Pp/W is the power available to the propulsion system. It can be shown 
that net thrust for any number of propellers can be determined with this method 
if: Ap is total propeller-disk area; Pp is total power absorbed by the pro- 

pellers; and all propellers have the same value of both tip speed and power 
loading (Pp/Ap). 

The program of Appendix A imposes both an upper and lower limit on the 
power-to-wei ght ratio of equation (18). As shown in the sketch below all 
energy not used for propulsion is stored. 
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Winds 

The model of the wind aloft is based on one set of wind data. This data is 
the 99 percent profile of reference 5 in figure 3 which describes a wind pro- 
file that is exceeded only one percent of the time at five sites in the United 
States. Figure 3 shows that the profile shifts substantially when the proba- 
bility of including all winds is decreased from 99 to 95 percent. The second 
99 percent profile of figure 3 is based on data from a world-wide set of sites 
(ref. 16) and indicates that the reference wind profile is reasonable but 
generally conservative. In the computer program of Appendix A, the magnitude of 
the wind at a given altitude is the product of the associated value from the 
reference profile and an amplitude factor . In the program, the direction 

of the wind vector at any altitude is assumed to remain constant at the same 
azimuth throughout an entire climb-glide cycle. 


DISCUSSION OF RESULTS 

The results of calculations of HAAP system performance are presented to 
show the effect of variations in aerodynamic, power-system, and other para- 
meters. Although not all of the combinations of values represent reasonable 
systems, the more extreme sets help to define trends. In most cases, the 
results are compared with the performance of the baseline HAAP system described 
in Table I. (This baseline system is similar to, but not indent! cal with, that 
of reference 9.) 

The variety of potential uses for a HAAP system has led to the use of 
several measures of performan'e to define results for most sections of this 
study. Requirements for a no;riinal ground-track pattern and the availability of 
sites for ground stations could produce emphasis on long endurance (total time 
per flight cycle) and zero-wind range. The need for high resolution in 
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observation tasks and wide-area coverage in communications tasks may produce 
some differences in the specifications for maximum altitude capability. Wide 
variations may also occur in the level of stored energy required to operate 
each payload as well as guidance and control systems. Therefore, the results 
presented for each parametric variation usually include range, endurance, final 
climb altitude, and stored energy. 


Typical Flights 

One cycle of flight is presented in figures 4 and 5 for each of three HAAP 
configurations with differing wing loadings. The essentially linear flight 
profile for gliding flight in figure 4 is a direct function of L/D. As shown 
in figure 5, the climb segment is affected by numerous parameters. The low 
wing-loading case of W/S = 50 Pa (1.0 lbf/ft2) has a fairly simple calculated 
history. During climb, the flight path is fairly linear, the rate of climb 
nearly constant, and the propeller provides thrust all of the time. The 
highest wing-loading case of W/S = 250 Pa (5.2 lbf/ft2) has an "s-shaped" 
climb profile (fig. 4) and climb history (fig. 5). Initially, the relatively 
smaller wing for W/S = 250 Pa does not receive enough power to start the 
motors. The airplane continues to store the received energy and to lose alti- 
tude as it glides nearer to the ground station (fig. 5(b)). When it is close 
enough to receive adequate power, the airplane uses all available power to 
climb. Near the end of the nominal climb period, the power received again 
drops below the minimum level. The airplane then glides and stores the 
received energy again. This latter glide segment illustrates the reasons that 
he can be less than the maximum altitude achieved in climb. 

Figure 5 indicates that there are strong relationships between performance, 
power available, and the flight path (defined with respect to the ground 
station). Consequently, the flight profile could be changed to maximize stored 
energy or some other parameter. All subsequent sets of results, however, are 
obtained for the simple, nominal type of flight at constant lift coefficient. 


Equivalent Airspeed 

The design value of Vg is important for a HAAP vehicle. As shown in 
equation (7) and figure 6, Vg is a function of both Cl and W/S. Maximizing 
L/D to improve range leads to the selection of the highest value of Cl that 
allows some margin of safety against stall. Requirements for adequate rectenna 
area and for long endurance (i.e., slow descent rate) can produce design anpha- 

sis on low values of W/S. Figure 6 shows that these design trends lead to low 
values of Vg. 

The effect of winds produce constraints on the minimum acceptable level of 
Vg. Figure 7 presents a wind profile that is exceeded only one-percent of the 
time at five sites (ref. 5), As indicated in a subsequent discussion on wind 
effects, this profile can provide a reasonable design criteria for HAAP vehi- 
cles that must avoid being blown away from a given site. The data suggests 
that Vg above 16.6 m/s (32.2 knots) is required if flight profiles extend to 
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as low as 18 kilometers in altitude. Application of this criteria to the data 
of figure 6 limits Cl as a function of W/S. 


Airplane Aerodynamic Characteristics 

The effect of W/S and L/D on HAAP performance are evaluated in 
figure 8. The parameter W/S also affects the power and propulsion system, 
since the rectenna is assumed to cover all the wing area S. Thus, decreasing 
W/S increases available power per unit weight. Large propulsion systems can 
then operate the propeller continuously at full power during climb. As shown 
in figure 5, this effect can result in a sustained high rate of climb. 

Figure 8 shows that reductions in W/S produce substantial improvements in 
attainable altitude and, below about W/S = 100 Pa, large increases in stored 
energy. Variations in L/D have relatively less effect on altitude and energy 
performance than variations in W/S. Range and endurance are both increased by 
reduced W/S or by increased L/D. (In the case of the baseline HAAP, Reynolds 
number effects change L/D as a function of altitude; however, the performance 
in that case can be shown to vary less than one percent from the performance 
for L/D = 45 and W/S = 100 Pa (2.1 Ibf/ft'^).) 

The effect of L/D can also be considered in light of the independent 
effects of Cl and Cq (fig. 9)* The value of Cg is assumed to be 
calculated as 

^D " ^D,o CL^/ftAe 

The set of Cg^o values used provides reasonable agreement between the maximum 
values of calculated L/D and those obtained from references 17 and 18. The 
effect of induced drag reduces the level of L/D that results from simply 
increasing Cl (fig. 9(a)). Both Cl and Cg^o appear to have significant 
effects on range and endurance. All of the trends for changes aerodynamic are 
in agreement with those determined in reference 11. 


Gliding Flight 

A simplfied analysis of gliding-flight endurance can be accomplished with an 
approximate solution to the expression of reference 8 for glide time between 
specified altitudes. Appendix B presents the development of an expression for 
a glide-time parameter tgVe (D/L) which is independent of configuration. 

An approximation for air density as a function of altitude allows the 
glide-endurance equation to assume integrable form. Figures 10 and 11 provide 
a means of comparing that the appropriate function with the values given in 
reference 4. Figure 11 shows that the two density models are in good agreement 
between 18 and 23 kilometers (59,000 and 75,000 ft), which the is altitude 
range of interest. 

The glide-time parameter can be used to determine the relative endurance 
achieved by gliding between different sets of initial and final altitudes. 

Figure 12 presents the glide-time parameter as a function of initial altitude 
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and altitude decrements. The computed results show that a given decrement of 
altitude yields a larger glide time at lower altitudes. This occurs because 
the vehicle travels more slowly through the denser atmosphere at lower 
altitudes. 

Results of gliding-endurance calculations are compared in figures 13 and 14 
for the computer program of Appendix A and the closed-form solution of 
Appendix B. The computer program has the advantage of accounting for accelera- 
tion effects and of using a more detailed model of density variation with alti- 
tude. The figures show that agreement between the methods is best at low 
values of W/S, L/D, and he . If the acceleration correction factor is 
removed from the computer program, the computer yields glide times which are 
virtually identical to those given by the closed-form method of Appendix B. 


Power Transmission System 

Climb performance is strongly affected by numerous interrelated parameters 
that characterize the power transmission system. As shown in equation (15) 
these include P/S, R, and kp. Parametric variations are considered here even 
though further development of microwave technology may lead to revisions of 
equation (15). The review of the present results is simplified by presenting 
only climb performance since gliding flight has already been treated. 

The character of beam-range effects is controlled by the exponent kp in 
equation (15). As shown in figure 15, focused power is independent of range at 
kp = 1. For any value of kp > 0, the equation requires that received power 
increases indefinitely as r/R approaches 0. In the present study, the values 
of r/R range from about 0.4 to 1.0, and the effect of kp does not reach 
physically implausible proportions for kp = 1. In an real system the trans- 
mitter would have a finite value of beam intensity at zero range; beyond a 
given range some value of kp would model the beam attentuation. Thus, 
increases in only kp imply a disproportionately large increase in actual 
transmitter power. Due to the large value of R, r/R 4 1 during climb; since 
the effect of kp is amplification at those regions, power intensity is always 
equal to or greater than the reference value P/S. 

Climb performance is presented as a function of kp in figure 16 and 
several climb histories are presented in figure 17. Increasing kp appears to 
allow the value of final climb altitude to increase asymptotically to a maximum 
and stored energy to increase exponentially. Since r/R 4 1, increasing kp 
simply increases available power at the vehicle. The calculated results appear 
to be opposite to the effects that would be anticipated from an increasing 
decay of beam intensity with distance; however, the short ranges and the 
implied large increase in transmitted power overcome the effects of decay rate. 

Climb performance is also sensitive to reference range R and the power 
density at that range P/S. Increasing P/S leads to large increases in 
stored energy (fig. 18) and allows the vehicle to climb higher. However, as in 
the case of (Pp/W)max “ ^ W/N in figure 18, motor size can limit altitude, no 

matter how much power is receiveu. Similar trends are shown for increasing R 

in figure 19. Increases in R or P/S are also associated with large 

increases in transmitted power. 
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The initial range and altitude for beam intercept also affects climb 
performance. Figure 20 shows that beam interceptions at longer range permit 
higher altitudes to be attained. However, the trajectories of these higher 
flying vehicles can reduce the amount of stored energy per flight cycle due to 
the attenuation of received power with range. This attenuation and the 
decrease of density with altitude combine to determine vehicle ceiling. As 
shown in figure 21, both the rate of climb and the energy storage for the base- 
line configuration can be estimated to be negligible at an altitude of about 
29 km (95,000 ft.). 


Propulsion System 

The effects of variations in overall propulsion-system efficiency are shown 
in figure 22. The computer program of Appendix A determines n as a function 
of J and Cp from a conventional propeller-performance table (ref. 15). 

This tabulated data does not reflect any effects of high-altitude, low Reynolds 
number phenomena on propeller aerodynamics. This omission, and other simplifi- 
cations, may lead to optimistic predictions of propeller performance. The 
result of operating with degraded thrust Td is a nearly linear decrease in 
attainable altitude (fig. 23). This indicates that even a small degradation 
in propulsion-system efficiency translates into noticeable performance 
decreases. 

The effect of relative motor size is shown in figure 23. The parameter 
(Pp/W)max »'eflects not only the maximum power that the propulsion system can 
absorb, it also indicates the ratio of motor size to total vehicle weight. The 
largest value of (Pp/W)max considered here is twice that of the baseline con- 
figuration. The computed results show that increasing the relative size of the 
motor generally leads to decreases in stored energy and to increases in attain- 
able altitude until a maximum performance level is achieved. Beyond that point, 
increasing (Pp/W)max is detrimental to performance. This variation indicates 
that the optimization of propulsion parameters is a function of wing loading 
(and rectenna size). 

A review of the calculated flight histories leading to the results of 
figure 23 reveals that the variation in performance is related to both trajec- 
tory characteristics and limits on the minimum power required. The vehicle 
with the larger motor may have to glide closer to the ground station before 
receiving enough power to overcome starting loads and other constraints. The 
more powerful vehicle climbs faster and generally flys a higher trajectory as 
it passes over the ground station. The more powerful vehicle then reaches the 
minimum (Pn/W) condition and begins its glide phase sooner. Detailed design of 
a HAAP will apparently be sensitive to constraints on minimum and maximum 
motor power. 

The effect of two propeller parameters on climb performance is shown in 
figure 23 and 24. The baseline value of tip speed (172 m/s) appears to be a 
good selection (fig. 24), although performance appears to be fairly insensitive 
to small variations in that parameter until the tip speed encounteres compress- 
ibility effects. The area ratio S/An is a somewhat artificial parameter that 
is a convenient element in equation (i8). That measure of relative propeller 
size is also set at a good value in the baseline configuration (S/Ap = 2.65). 
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Winds Aloft 


Although winds aloft can greatly influence the success of any given 
mission, wind effects on HAAP design are difficult to quantify. The statisti- 
tn of basic wi od data (refs. 5, 16, and 19) must be properly evaluated 

to avoid developing excessively stringent design criteria. Wind profiles that 
are exceeded only one percent of the time probably provide adequate design 
guidelines. The winds that exceed those limits tend to be associated with 
large storms occurring at lower altitudes. These more detectible, lower- 
altitude phenomon may provide enough warning to make appropriate changes in the 
flight program, such as maintaining as much altitude as possible. In addition, 
the relationships of wind direction at different altitudes are not considered 
in most sources of data. Nonuniformity of wind direction at different altitu- 
des may make HAAP operations easier than predicted for uniform wind direction. 

Operational limits imposed by winds tend to affect HAAP operations at lower 

nlclnn^ffi’ ^ ''e - "’/s (19 kt), the selection of a 

design value of Vg for lower altitudes will ensure an adequate margin of true 

airspeed V at higher altitudes. Thus, operations need not be restricted to 
the nominal low-wind region of about 20 kilometers (66,000 ft). 


HMP operations with actual real-time wind profiles will be more complex 
than for the flight trajectories considered in this study. Profiles for mean 
wind values from reference 19 show consistent trends with altitude of different 
seasons in figure 26(a); however, the associated data of figure 26(b) show 
there is a considerable variation possible between the mean and instantaneous 
values. Below 18 kilometers altitude, the mean winds blow predominately from 
west ^0 ®3st, although the instantaneous value appears to vary considerably 
(fig. 26(b)). Data from references 16 and 19 clearly indicate that winds at 
la kilometers and above are typically much stronger in winter. Despite the 
evidence of complexity, this study models winds on the basis of the profile 
shown in figure 7 and on the assumption of uniform wind direction. The wind- 

affects only the magnitude of the nominal profile 
that ^ f‘l directly reflect the probability level of encountering 


Studies were conducted of the effect of wind-profile magnitude and wind 

ground track. The first cases to be considered 
are those for the baseline HAAP configuration with a headwind or tailwind over 

hSLwinJ"?!! (fig. 27). Increases in wind-profile magnitude for a 

headwind (P =0 ) reduce groundspeed and increase the amount of time spent in 
passing over the ground station. The additional energy available through the 
extended climb period produces substantial increases in attainable altitude; 

headwinds affect the glide for a longer period of time and can 
substantially reduce total range. The reverse relationships appear true for 

arlrkilomJSr^'^^Jh T ° at k^ = 0.97 because headwinds 

kilometers, the initial altitude, can become no stronger without blowing 
the vehicles away from the ground station. ^ 


As shown in figure 1(c), adjustments to vehicle heading can cause the vector 
summation of wind and airspeed velocities to produce the desired ground track 

14 



for the HAAP (for sufficiently low wind speeds). However, if the vehicle rec- 
tenna is not exactly aligned with the transmitting antenna, the use of linear 
polarization will result in a reduction in power-transmission efficiency 
(eq. (14)). The effect of parametric variations in wind conditions is shown in 
figure 28. As in figure 27, the absence of calculated results for a given con- 
dition indicates that the baseline HAAP configuration could not fly in those 
winds. Typical performance near to limiting conditions is shown in figure 28(a) 
for y = 45 and = 0.74. As winds approach limiting conditions, the vehicle 

spends a large part of its climb time in slowly making headway at the lowest 
altitudes (near 18 kilometers); power storage increases significantly, but final 
altitude decreases. Figure 28 shows that as amplitude of the wind profile 
increases, only tailwinds permit flight. In all cases, the unsuccessful attempts 
at flight were terminated by winds at 18 kilometers blowing the vehicle away 
from the ground station. 

Flight with more severe wind profiles would be possible for all wind direc- 
tions if the baseline configuration or flight plan were modified. Previously 
discussed results show that increasing the design value of equivalent airspeed 
could allow the vehicle to operate In the presence of stronger winds. Another 
solution would be to increase the value of minimum altitude. As shown in 
figure 7, the nominal wind profile for this study is more severe at the lower 
altitudes. Figure 26 shows that such data is representative. An alterante 
solution would be to accept the cost and complexity of circular polarization, at 
least for the transmitter. The relative benefits of the last two methods are 
suggested in figure 29. If stored energy is not a limiting factor, the restric- 
tion of the flight profile to higher altitudes appears to offer a simple, viable 
sol ution. 

Although turbulence and wind shear affect the development on HAAP design 
criteria, these effects are not considered herein. Some limited data on these 
phenomena at high altitude are availabe in references 20 and 21. 


CONCLUDING REMARKS 

A parametric study of performance has been conducted for remotely-piloted, 
microwave-powered, high-altitude airplane platforms. The nominal flight plan 
consists of climb and glide cycles: while receiving power, the vehicle climbs 

and stores excess energy; it then glides back down to a minimum altitude above a 
microwave ground station. 

Calculated results identifed several basic trends. Low values of wing 
loading and high values of lift coefficient were shown to result in long range, 
long endurance, and low equivalent airspeed. Wind effects constrain the lower 
limits of both equivalent airspeed and operating altitude. Calculations also 
showed that power-transmission and propulsion-system characteristics could 
strongly affect climb performance. An approximate, closed-form solution was 
developed to predict gliding endurance. 
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APPENDIX A 


COMPUTER PROGRAM FOR HAAP PERFORMANCE 


^ computer program has been developed to calculate the performance of a 
HAAP vehicle. This appendix contains a listing of the program, a sample input 
file and the corresponding sample set of output listing. The results presented 
in the output listing can be interpreted with the description of variable names 
given in Tables II and III. 


The program calculations and logic are based on the HAAP operating proce- 
dures as described in the main report. The program calculates the flight tra- 
jectory and system performance at specified intervals of time. These intervals 
are 10 seconds for climb and 20 seconds for glide unless the end of climb or 
glide is approached; at that point, the intervals are adjusted to be one-tenth 
their previous value. The only configuration change allowed during a given 
flight is the folding or unfolding of the propellers. 


The program yields results for parametric studies. The first set of output 
information describes initial conditions in terms of the characteristics of the 
airplane aerodynamics, propeller and power-system variables, wind and trajectory 
parameters. The listings presented in columns provide histories of performance 
and flight mechanics. For each run, the input parameter being varied is listed 
in the first column on the left. Each set of parametric variations may be con- 
ducted for performance at a single point (with respect to the ground station), 
during climb or glide only, or throughout an entire climb and glide cycle. 


The sample case included in this appendix illustrates the effect of wind 
magnitude. Performance is calculated for the baseline configuration HAAP with 
winds at nght angles to the nominal ground (y = 90°). The required inputs ar 
N1 - 3, N2 = 10, AI^U = 90., SI = 0., SF = 1.0, and SS = 0.2. Results indicate 

strength wind profile does not allow the vehicle to initiate climb 
at 18 kilometers. 
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i 


PROGRAM HAAP ( !N PUT* OUTP UT » TA P E ■ I N P UT» T A PE 6» OUTP U T) 


'D 


10 


15 


20 


25 


30 


35 


40 


r 

OIMENSIOM A(5)* B0N(17)* BDU7) ,Z(8>20) 

C A- ALPHANUMERIC LABEL* 3DN- NAMES OF ELEMENTS OF BASELINE DATA ARRAY 

C BD- BA5EINE DATA ARRAY* Z- FINAL OUTPUT ARRAY 

C 

COMMON /PAAH/ WOS*CL* 3LDD*HL0U* Ti*SOAP* POS*RR* PDWL*WK* AMU* 

1 XS*HS* PO^K*TDOT*HI#PKR» AK*ETA* GAHMA*POW* POWP*POWS> PSID* 

2 R*RLQD* ROC*THETA* VE*VG* VT*PCP* PJ*N4*N5*X*H 
EQUIVALENCE {BO(l)»WaS) 

C 

NAMELIST/DD/ WOS* C L* 3LQJ * HL QD * TS > SOAP* POS* RR * POWL * WK * AMU* XS* HS* 

1 3I»SF*SS*N1*N2*N3*N4*TD0T* POWR*HI*RKR 

C 

DATA WOS * C L* BL3D *HLOD*TS* SOAP* POS* RR* POWL *WK* AMU*XS*HS* 

1 SI*SF*SS*N1,N2*N3*N4*TDCT* POWR*HI*RKR/ 

2 144. * 0.9* 36. 6*. 41 8* 172. *2. 653 *1.1* 50. *8. 6 2* 0. 0*0.0* 

3 40. * 18. 0* 0.*0 .*0. *2*2*5 0* 1* 1.0* 0.25*2 5. *1./ 

C 

DATA BDN/7H W/S* 6H CL* 8H L/D-B* 8H L/D“H*6H TS* 

1 8H S/A-P* 7H P/S* 6H RR* 8H P/W-L* 6H WK* 

2 6H MU* 5H XS* 5H HS* 9H P/W M-H* 7H TDOT* 

3 6H HI* 6H PKR/ 

C 


C 

r» 

BASELINE DATA ARRAY 

1. W/ S 

2. CL 

3. 

L/D-B 

4.L/D 

-H 5.TS 

6. S/A-P 

C 

C 

C 


7. P/ S 

3.RR 

9. 

P/W-L 

10. WK 

11 .MU 

12. XS 


13 .HS 

14. 

POWR 

15. 

TDOT 

16. HI 

17. RKR 


C 

C parameter variation code si- initial VALUE 

C Sr- FINAL VALUE 

C SS- STEP INCREMENT (POSITIVE OR NEGATIVE) 

C 

C CONTROL CODE N1-(5INGLE POINT* CLIMB* TOTAL FLIGHT* GLIDE ONLY) 

C N2-(ELEMENT IN ARRAY BD TO BE VARIED) 

C N3-(NUMBEn of CALCULATION POINTS PER LISTED LINE) 

C N4-(RECTENNA POLARIZATION- LINEAR OR CIRCULAR) 

C N5-{C0DP=1 WHEN VG< 0) 

C 

100 FORMAT { 1H1»5X*5A10// 5X , 14HAIRCRAFT AERC.* 7 X » 9 HP RO PE L L E R * ilX* 

1 5HP0WER* 14X*5HWINDS* 9X*11H6TART POINT* 4 X * 12H V AR I A BLE SET* 



00 


2 9X,4HCHDP// 5X,4HW/S»>F6.1, 5H N/!i2^ 6X, 3HTS»,*F6.i# 4H H/S>7X, 

3 4H?/S*»F5,2> 6H KW/M2> 4X> 3HWK«»F5.2^ 6X# 3HXS-»F6.2^ 3H KM# 4X# 

45 4 6HFIR3T = » PB.3# 5X# i 1 4X # 3HN 2» # 1 3 / 6 X # 3HC L =# F 5 . 2# 9X# 

5 6HS/A-P«#F6 .3# 12X » 3HRK=» F 5 , i# 3H KM# 7X# 3HMU«#F5.1# 4H DEG# 2X# 

6 3HHS«#F6.2# 3H KM# 4X# 6HF I.MAL »# F8. 3# 5X# 3HN3«»13/ 114X# 3HN4-# 

7 13/ 5X»4HL/D = # F5.1# 2RX# SHMAX R/w=#F5.2# 6H KW/KN# 18X# 3HHI-# 

e F6.2# 3H KM, SHST'^P*# F8.3# 6X# // 2X# 7HL/D(H)-, 

50 9 F6.3#27X,BHKIN P/W»# F4.2# lOH X MAX P /W ///) 

101 F3RMAT (5ftlO) 

111 FTP ''AT ( 2X# A10#4X#1HX#5X#1HH# 5X# 3HR / C # 3 X # 5 HP / W-F # 2X# 5HP/W-S#2X# 

1 dHGAMMA# 2X» 5riTHSTA#4X#lHR#6X,2HVG# 4 X# 2 H VT #5 X # 3H V E C # 5X # 

2 1HT#6X# 2HAK#4X# 3HETA# 3X# 2HCP# 4X# IHJ #5X# 3HPSI// 

55 3 16X#2HKM# 5 X# 2H< '1 # 4X » 3H M/ 3 # 4X#3HW/N#4X» 3HW/N#4X# 

4 3H0EG# 4X»3H0EG# 5 X #2HKM# 4X# 3H M /S # 4X#3HM/S# 4X#3HM/S# 

5 4X,3HSEC»29X# 3HDFG // ) 

113 FORMAT (2X#F9.5#1X#»6F7.2»F7.1#F7.2#3F7.1#F3.0#2F6.3#2F6,3#F6,1) 

114 FORMAT (5X/) 

oO 115 FORMAT (2X,F9.5#1X# 3F7.2# F21.2# 14X# 3F7.1# F8.0# F8 .4 ) 

200 FORMAT (//2X»A10# 4X,2HXC» 5X#2HHC# 4X#2HTC# 4X#5HE/W-S# 3X# 

1 5HE/W-T, 4X,2HXT. 7X#2HTT// 16X # 2HKM# 5X#?HKM# 4X#2HHR# 

2 4X#4HKJ/N, 4X.4HKJ/M# 5X#2HKM# 7X#2HHR/) 

201 FORMAT (110X#3HTG«# F6 . 3 # 3H HR //) 

b5 202 FORMAT ( 3X , F 8 . 3 # 1 X# 2F7 . 2 # F 6 . 3 # 2 F 6 . 3 # 2 F 8 . 2 ) 

r 

c 

GO TO 3 

2 X» 0.0 

70 H« 

IF (Ml .EQ.1.0R.N1.EQ.4) GO TO 3 
WRITE (6# 200) 3DM(N2) 

DO 300 I*1#M 

WRITE (b# 2 02) Z ( 1# I ) # Z ( 2 # i ) # Z ( 3 # 1 ) # Z { 4# I ) # Z ( 5 » I ) # Z { 6# I ) # Z ( 7# I ) , 

75 1 Z(8#I) 

300 CONTINUE 

3 READ (5# mi) A 

IF { EOF (5 ) ) 1 09, 4 

4 READ (5#DD) 

Uo W-ITF (6»100) A,VDS»TS#POS#rtK,Xj#Sl»Ml# N2# C L # S □ A P# RR# AMU # HS # SF # 

2 M3#N4# BLJO# fowl# HI#3S# HLDD#P0WR 

0 

C N:2- f'ARAMFTFR FOR VARIATION# CHOSEN FROM ARRAY BO 



83 


90 


93 


IOC 


105 


110 


115 


120 


C 


C 


c 

c 


c 

c 


c 

c 


7 WRITE (&»111) BDN(N2) 

INITIALIZE PARAMETER NEW STARTING POINT IS NEW X OR Y 

N3= 0 

BD(N2)» SI 
GO TO 17 

13 Z (2,N)« X 
Z (3^N )«= H 
Z (A,N) = TT 
Z ( 5» N >1= EOWS 
Z (6>N)« EOWT 
IF {N1.EQ.3) GO TO AO 

INCREMENT PARAMETER 
lA CONTINUE 
N5 ■ 0 

I? (N1.EQ.3) Z(7.N)» X 
IF (N1.E0,3) Z( BjN) » TG 
IF (Nl.EQ.A) WRITE (6»201) TG 


16 BD(N2)= 3D(M2) +SS 
DELTA* SF -BD(N2) 

IF (SS.GT.O. .AND. DELTA. GE.O.) GO TO 17 
IF (S5. LT. 0. .AND .DELTA. Lt . 0. ) GO TO 17 
GO TO 2 

17 N* N+1 
Z(1»N)- BD(N2) 

GO TO (20» 30, 30,39) ,H1 

CALCULATE VALUES AT ONE ROINT 
20 T= 0. 

X* c. 

H» HS 
EOWS= 0. 

EOWT= D. 

GAMMA* 0. 

CALL RCLIMB 

WRITE (6,113) BD(N2), X,H, ROC,POwP, ROWS, GAMMA, THETA, R, VG,VT, 
1 VE»T, AK,ETA, ?CP,PJ, PS ID 

WRITE (6,11A) 

GO TO lA 


125 



130 


135 


140 


C 

C CALCULATF TOTAL CLIMB PHASF 

30 NK» 0 

PQW1= 0.0 

P0WS1» 0. 

GAMMA- 0. 

THFTA* 0.0 
T» C. 

EOwS- 0. 

EOWT« 0. 

RR95» RR* .95 
TT- 0. 

X- 0.0 

TDEL« 10. 

WRITE (6,114) 


145 


150 


155 


160 


165 


C N3 - PRINTOUT INCREMENT FOR COMPLETE CLIMB 

31 IF (NK.EQ.N3) NK« 0 
GAMMA- GAMMA/57. 2957 

32 CALL RCLIMR 

IF (N5.GE.1) GO TO 14 

IF (P .GT .RR95. AND. THETA. LT .50. ) TDEL- 1. 

IF (T.EQ.O.) GO TO 35 
NK- NK+1 
C 

C X AND H GIVEN IN KM 

33 X* X +VG+TDEL/1000. 

H- H +K0C*TDEL/1000. 

C 

C SPECIFIC-ENERGY INCREMENTS FROM AVERAGED POWER FOR TIME INCREMENT 

C STORED SPECIFIC ENERGY (E/W-S) AND TOTAL UTILIZED SPECIFIC 

C ENFRGY (F./W-T) ARF GIVEN IN KJ/N 

EOWS- EOWS +(P0WS+P0WSl)*TDEL/2000. 

POWSl- POWS 

EOWT- EOWT +(POaI + POWl) *TD EL/^000 . 

POWl- PCW 

IF (R . GT.RR.AND.THETA.lt. 90) GO TO 35 
IF (X.LT.O.) GO TO 35 
IF (NK.EQ.N4) GO TO 35 
T» T+TDEL 
GO TO 31 



C WPITE DATA HNE INCREMEKT QF CLIHB OR FINAL CLIMB POINT 

170 35 CONTINUE 

WRITE (6>113) BD(N2)> RQC>POWP» P0W3>GAMMA» THETA#R, VG»VT# 

1 VE»Tf AK»ETA, PCP>PJ» PSID 

TT= T/3600. 

T= T+TDEL 

175 IF (X.LT.O.) GO TO 13 

IF (F.GT.RR.AND.THETA.Lr.90) GC TO 13 
GO TO 31 

r» 

c 

18u C CALCULATE GLIDE PHASE 

C 

39 T» 20. 

TOEL= 20. 

X= 0. 

185 H» HI 

AG NG = 0 

HSIOI- 1.01* HS 

T- T -TDEL 
PDF= 2.653/SDAP 

190 C PDF IS PROP DRAG FACTOR-PROPORTIONAL TO RATIO OF DISK AREA TO WING AREA 

TDEL « 20. 

GAMMA- 0.0 

A1 GAMMA- GAMMA/57.2957 
C 

19t; C BEGIN CALCULATION FOR NEW GAMMA AT NEW ALTITUDE 

A2 PLOD- BLGD +HLOD*H 
IF (HG.EQ.N3) NG* 0 

C DECREMENT L/D DUE TO DRAG OF FOLDED PROPELLER 

IF (Nl.LT.A) RLOD- RLOD -1.5+PDF 
200 KK- 0 

r> 

U 

C ITERATE FOR GAMMA 

A3 KK= KK +1 

IF (KK.EQ. 10) GO TO 50 

205 VE- 1.27775*SQRT ( WOS* C 0 S ( G aMM A ) / C L ) 

CALL ALTF (AMU, VE, WK, H, PSI, VT, VG, N5) 

VY= -VT*C OS (GAMMA) /RLOD 

CALL ACCSL(VY,H,VG, A MJ , V E , WK , G AMM A, A K) 

ROC- VY/ ( 1. +AK ) 

210 G CALCULATE RESULTING CLIMB ANGLE 


ro 



215 


220 


225 


230 


GAMMAC- ASIN(R0C/VT) 

DELG* AB5 ( GAMMAC-GAHMA) 

IF (DELG.LT ..0002) GO TO 5C 
AJUST CLIMB ANGLE AND REPEAT 
GAMMA* GAMMAC 
GO TO 43 
50 NG * NG + 1 

X AND H GIVEN IN KM 

54 X* X +VG’!'TDEL/1000. 

H» H +ROC*TDFL/1000. 

T» T + TDEL 

IF (H.LT.HSIOI) TOEL« 2. 

IF (HG.E0.N3) GO TO 55 
IF (H.GT.HS) GO TO 42 

55 GAMMA- GAMMA+57. 2957 

WRITE DATA FOR ONE INCREMENT OF GLIDE OR FINAL POINT 
WRITE (6,115) B0(N2),X,H, ROC, GAMMA, VG, VT, VE, T, 
IF (H.GT.HS) GD TO 41 
TG= T/3600. 

Z(7,N)« X 
Z(8,N )» TG 
GD TO 14 


235 


C 


199 STOP 
END 



1 c 

c 

SUBROUTINE DENSITY (CH» SIGNA) 

C 

p C CURVE PIT TO 62 ATMOS. F GK CALCULATION GF DENSITY RATIO 

C 

C INPUT: ALTITUDE IN KM; OUTPUT: OIMENSIONESS DENSITY RATIO (SIGMA) 

DIMENSION DCl(15),DC2(lP) 

C 

10 C sigma- E **(CC1«H +CC2*H«*2) WHERE H IS IN KM 

C 


DATA (DCK I) »I-1#1 5) / . 0956 5 54# .0948554, .0955 529, .0960089, .0942258, DCl 

1 .0942258, .0770834, .0879373, .0962238, .1027082, DCl 

2 .1045655, .1107329, .1160101, .1204581, .1243220/ DCl 

15 DATA (DC2(J),J«1, 15)/. 117337, .117337, .124898, .133965, .143754, DC2 

1 .143754,. 307572, .230044, .178253, .142229, DC2 

2 . 132942, .104908, .082920, .065612, .052013/ DC2 


C 

ICH* 1+ IEIX(CH/2. ) 

20 C 

CCl- DCKICH) 

CC2- DC2{JCH) 

C 

IF (CH.LE. 11. . 0R.CH.GE.12. ) GO TO 20 
25 CCl- .0675418 

CC2= .387035 
C 

20 SIGMA- EXP (-CC1=1‘CH -CC 2^C H*CH/100. ) 
RETURN 

30 END 


ro 

CO 



ro 

45 » 


SUBROUTIolt ALTF( AMU»V£# WK»H»PSI > yJrMGfHb) 

C CALCULATE TRUE AIRSPEED# wINDSPEED# AND GROUNDSPEED - SI UNITS 

5 C inputs: wind AZIMUTH# EQUIVALENT AIRSPEED# WIND SCALE FACTOR# AND 

C ALTITUDE; OUTPUT: GROUND-TRACK OFFSET ANGLE# TRUE AIRSPEED# AND 

C GROUND SPEED. (ALL SPEEDS IN M/Sj ALL ANGLES IN DEGREES.) 

C FOR WK= 1.# resulting wIND PROFILE IS FOR 99K INCLUSIVE 

C PROFILE FOP 5 LAUNCH SITES FROM NASA TM 78118. 

10 C 

300 FORMAT (AX# 2AHWIND SPEED TOO LARGE AT #FA.l# AH KM.) 

CALL DENSITY (H.SIGHA) 

VT« VE*( SIGMA)«*(-.5 ) 

11 IF (H.GE. lA. ) GO TO 50 

VWe WK^88. 

GO TO 62 

50 IF (H.6F.15.) GO TO 51 

VW= wKt (8P. -18.*(H -lA.)) 

20 GO TO 62 

51 IF (H.GE. 20.) GO TO 52 

VW« WK+ (70. -5.8*(H -15.)) 

GO TO 62 

32 IF (H.GE. 23.) GO TO 53 

25 VW« WK+Al. 

GO TO 62 

53 VW« WK’>= (Al. +A.7778+(H -23.)) 

62 SPSI= VW*SIND( AHU)/VT 
IF (SPSI.GE .1. ) GO TO 6A 
30 PSl* ASIN(SPSI) 

VG* VT*C3S(PSI) -VW*COSD (AMU) 

GO TO 65 

6A WRITF (6#300) H 
N5= 1 

35 05 CONTINUE 

RETURN 

END 



1 


t; 


1 r\ 


ib 


20 


25 


30 


35 


40 


U 

r- 

L 


U 

U 

c 

c 

c 


5UBPGUTIME PPOPCAL (PCP>PJ^ETA) 

DT^^tNSIDN PT(15,41) 

INPUT: P-^OPFLLEP POWEP COEFFICIENT AND ADVANCE RATIO 
output: propeller EFFICIENCY FACTOR 

data SOURCE: HAM. STD. CHARTS FOR CL-I* 0.3» AF« 80. » AND THREE BLADES 
EACH DATA STATEMENT GIVES VALUES OF ETA AS CP* RANGES FROM .0 TO .35 


data {PT(I» 1)>I = 1*15)/.0* .74>.66*.57*.48».42*.37*.33».28*.2^*.22> 
i.l9» . 17* . 16* .14/ 

CAT A (®T(I* 2)*T«1*15)/. 0* .78* .71* .61*. 54*. 47*. 41* .36* .32* .20*. 2 4* 
1.22* .19*. 175* .16/ 

DATA (PT(I* 3)*I*l*15)/.0*.80*.7 6* .66* .5 6* .52*. 45* .41*. 35*. 31 *.2 7* 
1.24*. 22*. 19*. 18/ 

DATA (PT(I* 4)*I»l»15)/.0*.83». 80* . 7 1* . 6 2* . 5 6* . 49 * .45* . 39 * . 35* .3 0* 
1.27* .24* .22* .19/ 

DATA (PT(I* 5)*T“l*15)/.0*.84*.82* .7 5*. 67*. 60 *.54*. 48*. 43* .38* .33* 
1.29* .26*.24*.22/ 

DATA (?T(I» 6)*I «1 *15)/. 0*. 85*. 84*. 78*. 70*. 63*. 57*. 52*. 46*. 41*. 37* 

1 .32*. 28*. 25, .24/ 

DATA (PT(I* 7)*I»l*15)/.0*.87*.86*.80* .73* . 67* .61, .55* .50* .44* .40* 
1 .35*. 3 1*. 23*. 25/ 

data {PT(I* 8)*I»l* 15)/.0* .88*. 37 *.82*. 76 *.70*. 64* .59*. 53*. 47*. 4 3* 

1.38. . 34*. 30* .29/ 

DAT A (PT(I* 9)*I»l*15)/. 0*. 88*.88*.84*.78*.73*.67* .62* .57, .52* .46* 
1.40* .37*. 34*. 29/ 

DATA (PT(l*10)*I = l*15)/. 0*.88*.9 9* .8 6*. 81*. 75* .70* .65* .60* .55* .49* 
1 .44*«39*«36* .33/ 

DATA (PT(I*ll)*I»l*15)/.0* .88*. 90* .8 7*. 8 2*. 78*. 72*. 68*. 63*. 58* .52* 

1.47. . 42* .39* .34/ 

DATA (3T{I*12)*I*‘1*15)/.0*.87*.9 0*.8 8*.84*.80*.75*.70*.65*.61*.56* 
1 .5* .45* .41 * .38/ 

DATA (PT{I*13)*I*=i,15)/.0*. 67* .9 14* .89* .85*. 81 *.77, .72*. 68*. 63*. 6* 
i. 54*. 4 8 *.45* .40/ 

DATA (»T(I*14) *1*1*15)/. 0* . 80 * .913*. 9*. 87*. 83* .73*. 74*. 7*. 6 6*. 62 * 
1 .56* . 51* . 47* . 44/ 

DATA {?T(I*15)*I*l*15)/.0,.35*.915*.901*.88*.e4*.b*.76,.72*.68* 

1 . 64 * . 6 0, . 5 5* . 50* . 47 5 / 

DAT A {PT(I*16)*I*1*15)/.0*.64*.916*.91*.89*.85*.32*.78*.74*.7*.67* 
1 .62* . 575, . 54* .49/ 


.40 
.40 
.45 
.45 
.50 
.50 
.55 
.55 
.60 
.60 
.65 
.65 
.70 
.70 
.75 
.75 
.80 
. 80 
.85 
.85 
.90 
.90 
.95 
.95 
1.00 
l.CO 
1.05 
1.05 
1.10 
1.10 
1. 15 
1. 15 


ro 

cn 



ro 

CT> 


H 


i>Q 


55 


60 


65 


70 


75 


80 


DATA (PT (I, 17)> I«l, 15)/ ,0> .84, .917, . 912> .8 95,. 86, . 83>. 79, . 7 5, .73, 

1 .69, .65, .60, .55, .53/ 

data (PT( I, 18), 1»1, 15)/. 0,. 83,. 918, .915, .9, . 87,. 8 4,. 81,. 77, .74, .7, 

3«67>.63>.58>.55/ 

DATA (PT(I,19),I«1,15)/.0, . 82,.917,. 916, .902, .88, .85, . 82, . 79, . 75, 

1 • 73 ^ 9 69 f 0 66f m bl f .SB/ 

DATA ( PT (I, 20), I»l, 15)/. 0,. fa, .916, .92, .91, .89, .86, .83, .8, .77, .74 , 
1 •?!> • 675> • 64^ • 60/ 

DATA (PT (1,21 ), I »1, 15) /.O, .a, .915, .92, .912, .89 5, .8 7, .8 4, .81,. 78, 

1 • 75 f • 7 2 p • 7 Gp 0 67 p *63/ 

DATA (PT(I,22),I-1,15)/. 0,.8,.915, .92,.914,. 9,.33, .85, .82, . 8,.77 , 

l*7^^*71^#6Bf #65/ 

DATA (PT(I,23),Ial,15)/.0,.8,,9i5, .9 2, .916, .90 5, .8 85, .86, .83, .81, 

1 •78^ %7 5 9 0 73 p *70# *67/ 

DATA (PT(I,24),I»l,15)/.0,.8,.9i2,‘. 9 2, .918,.91,.89,.8 7, .84,.82,.8, 
1*77. ♦ %7 ^ p 07 2 p 0 69 f 

data (PT( I, 2 5), T» 1,15)/. 0, .6, .91, .92, .92, .911,. 90, .88, .85, .83, .6 , 
3, «76> #755^ «73> • 70/ 

DATA (pT(I,26),I»1,15)/.0,.8,.91,.9 2, .9 2, .914,.9, .885, .86, .8 4, .81, 
1 .79,. 77, . 74, .72/ 

DATA ( PT (1,27) ,I»1,15)/ . 0, .79, .905, . 92, .92, .915, . 905, . 89, . 8 7, .84 5, 
1 .82, .80, .78, .755, .735/ 

DATA (?T(I,28), I«l,15 )/.0,.78,.902,.92, .92, . 917, . 909, . 89 5, .8 75, 

1 .85, .83, . 81, .79, .77, .75/ 

DATA ( PT ( 1,2 9), I« 1,15)/. 0, .77, .9, .92, .92,. 9 17, .91, .90, .8 8, .86, .84, 

1«82>»60>»7B^»76/ 

DATA (PT (I ,30 ),I»1, 15) / . 0, .77, .9, . 917, . 9 2, .9 2, .913,.903, . 89, .8 7, 

1 .85, .83, .31, .79, .77/ 

DATA ( PT( 1,31 ) ,I«l,15)/.0, .76, .89, .915, .92,.9 2, .915, .907, . 893, .88, 
1 .86, , 84 , . 82, . 80, .78/ 

DATA (‘’T(l,3?),I*l,15)/. 0,. 75,.88, .914, .92,. 92,.-916, .909, . 895, . 88, 

1 *36>*85^#83>#81f*79/ 

DATA (PT(I,33), I-l,15)/.0, .74, .87, .911, . 9 2, . 9 2 , . 91 7, . 912, . 9, .89, 

1 0^7 p 0^55 9 0 635 p 0 32 p 0 b 0/ 

DATA (PT( I, 34), 1 = 1,15) /.O, .72, .86, .91, .92, .92, .9 19,. 914, .903, .89, 

1 •87^«B6>*8^5>«325>i31/ 

DATA (PT(I, 35), 1=1, 15)/. 0,. 7, .86, .90 5, .919, . 92, .92, .915, .90 8, .9, 

1 •8B>»‘3 6 5^ •BSf »B3^ *82/ 

data (PT(1,36),I»1,15)/. 0, .7, .85,. 9, .917, .92, .92, .917,.91,.9,.89 , 
l*87>i66^ *8^^ 0 62 5 f 

data (PT ( T,37) , 1 = 1, 15) /. 0, .7, . 85, .9, .915, .92, . 92, .918, .912, .902, 

1 .89, . 88, . 865, . « 5, .83/ 


1.20 
1.20 
1.25 
1.2 5 
1.30 
1. 30 
1.35 
1.35 
1.40 
1.40 
1.45 
1.45 
1.50 
1.50 
1.55 
1.55 
1.60 
1.60 
1.65 
1.65 
1.70 
1.70 
1 .75 
1.75 
1.80 
1.80 
1.85 
1.85 
1.90 
1.90 
1.95 
1.95 
2.00 
2.00 
2.05 
2.05 
2.10 
2.10 
2.15 
2.15 
2.20 
2.20 



85 


90 


D4TA (0T(I#36)^ 1 = 1, 15)/.0> .7».85#.89,.913».92>.92#.913f.9l3».907» 2.25 

1 .695, .835, .87, .855, .84/ 2.25 

DATA (dT(T,39),i«i,15)/.o, . 7, .84, .88,.911,.92, .92, .91 8, . 91 5 ,. 908, 2.30 

1 .9 , .89,. 875, .36, .845/ 2. 30 

DATA (PT(I,40),I=1, 15)/.0,.7, .83, .86, .91, .918, .92,.918,.9i5,.9i, 2.35 

1 .902, .89, .88, .865, .35/ 2.35 

data ( PT ( 1, 41) , I»1, 15) / .0, .7, . 82, . 8 7, .9 07, .916, .918, . 918, .916, 2. 40 

1 .912, .905, .895, .885, .87, .86/ 2.40 


95 


100 


105 


C INTERPOLATION OP ETA FROfI HAM. -STD TABLES FOR AF» 8 0, B»3 , CL I « . 3 

RCP= 40.*PCP +1. 

ICP» IFIXCPCP) 

DCP= PCP -FLOAT( ICP ) 

RJ« 20.*PJ -7. 

IJ« IFIX(RJ) 

DIJ= RJ -FLOAT{ IJ) 

C POINTS A & Q AT GIVEN C? VALUE; POINT A AT LOWER J VALUE THAN POINT B 

PTA» (1.-DCP)^PT(ICP,IJ) + PT ( 1C P+1, I J) P 
PTB« (l.-DCP)^PT(ICP,IJ+l) +?T( 1CP + 1,IJ+1)^DCP 
ETA* PTA +(PTB-PTA)-»'OIJ 
RETURN 
END 



ro 

CO 


1 


5 


10 


15 


20 


C 

C 

SUBROUTINE ACCEL ( V Y» VG> AMU» VE # WK » GAMM A> AK> N5 ) 

C INPUT: VERTICAL VELOCITY IN M/S, ALTITUDE IN KM, GROUND SPEED 

C IN M/S, WIND AZIMUTH IN DEG, EQUIVALENT AIRSPEED IN M/S, WIND SCALE 

C FACTOR, AND FLIGHT PATH ANGLE IN DEG; OUTPUT: ACCELERATION CORRECTION 

C FACTOR 

IF (VY.LT.O) GO TO 8A 
Yle H+.l 
Y2» H-.l 
GO TP 85 
bA Yl= H-.l 
Y?* H+.l 

85 VAV« SORT (VY*VY +VG=«'VG) 

CALL ALTF ( AMU , V F, WK , Y1 , PS 1 1, VT 1, VGl, N 5 ) 

IF (N5. EQ.l ) GO TO 87 

CALL ALTF (AMLI,VE,WK,Y2,PSI2,VT2,VG2,N5) 

IF (N5.E0. 1) GO TO 37 

Vl« SORT (VG1**2 + (VT1*S IN (GAMMA) ) +*2) 

V2» SORT (VG2**2 +{ VT2*S IN (GAMMA) ) »t2) 

DF.LV= V1-V2 

AK* DELV=!‘VAV/1950. 

87 CONTINUE 
RETURN 

end 


25 



1 


c 




10 


15 


20 


25 


30 


35 


40 


SUBROUTINE RCLIMB 

COMMON /PA4H/ kOS»CL» BLQD»HLOD, TS^SOAP# POS»RR, POWL^WK^ AMU# 

1 XS#HS» P0WR#TDL1T, Hi#PKR, AK,ETA/ GAMMA# POw# PQWP#POWS# PSID# 

2 R»PLOD# ROC#THETA, VE#VG# \/T#PCP» PJ#N4,N5#X»H 
DATA Cl/ 9.9397/ 

C Cl EQUALS (PI=«=*4 )/(8.(S.L. DENSITT)) 

r 

C BASIC PARAMETERS 

XR= X -XS 

R» SORT (XR«XR +H«H) 

THETA* ATAN2(H,XR) 

PLOD* BLGD +HLOO^H 
KODE= 1 
KK* 0.0 

60 VF* ] .27775«SQRT{W0S*C03{GAMMA)/CL ) 

C NOTE: VE IS CORRECTED FOR FLIGHT PATH ANGLE# GAMMA 

CALL ALTF ( AHU#\/F#WK#H,PSI#VT# VG» N5) 
ic (N5.EQ.1) GO TO 90 
CALL DENSITY (H#SIGMA) 

C 

C CALCULATION OF POWER - R E C E I V E D# AV A1 LAB L E AND STORED 

C 740 FACTOR IS 1000 W/<W X .74 EFFICIENCY FACTOR 

ANGLE* 3.1415926 -THETA +GAMMA 
7 0 POw=( {RR/PO-X+R KR)* 740 . ( P 0 S / WO S ) * SI N ( A NGL E ) 

IP (N4.EQ.1) PDW* PQW« ( C 0 S ( P S I ) ) «*2 
KK* KK+1 

IF (KK.GT.IO) GO TO 90 
POWERL= ODW/PDWL 
IF (POWFRL.GT. POWR) GO TO 75 
C KEED PROP FOLDED AND STORE ALL INCOMING ENERGY 

FT A » 0.0 
POWP* 0.0 
p J = 0.0 
PCP= 0.0 

C DECREMENT L/D TO ACCOUNT FDR DRAG OF FOLDED PROPELLERS 

RLOD= RLDD -1.5 
KODE* -1 

p n w s * PON 

GO TO 33 

75 DPOW* PDW -PDWL 


l\5 



CO 

o 




50 


55 


60 


6 5 


70 


75 


80 


IF (DPQW) 76,76,77 

C ALL POWFi? TO PROP 

76 PQWP« ^DW 

POWF» 0.0 
GO TD 78 

C POWER TO PROP AND REMAINDER TO STORAGE 

77 POWP« POWL 

PQWS» OPOW 

76 IF (KOOE) 79,79,82 
79 RLOD» RLOD +1.5 

C 

C CALCULATION OF NOND I ME NS I ONAL CH AP AC TER ITl CS OF PROPELLER 

62 PJ« 3.14159*VT/TS 

POAP= POWP*WOS* SOAP 
PCP» C1*P0AP/(SIGMA*TS*^3) 

CALL PROPCAL (PCP,PJ,FTA) 

C 

C CALCULATION OF RATE OF CLIMB - THRUST AND DRAG COMPONENTS 

C TDOT IS RATIO OF ACTUAL, DEGRADED THRUST TD THRUST FROM TABLE LOOK-UP 

ETA= ETA^^'TOOT 
83 VYT=« ETA7PDWP 

VYO= VT*COS (GAMMA) /RLQD 
VY= VYT-VYD 

CALL ACCEL (VY,H, VG , AM U, V E , WK, GA MMA, AK) 

IF (N5.GE.5) GO TO 90 
RGC= VY/(1.+AK) 

C 

: CALCULATE RESULTING CLIMB ANGLE 

GAMMAC= AS INCROC /VT) 

DELG= A35 ( GAMMAC-GAMMA) 

IF (DElG.lt. .001) GO TO 90 

c 

C AJUST CLIMB ANGLE AND REPEAT 

GAMMA* GAMMAC 
GO TO 60 

C 

C CALCULATION FOP GAMMA (FLIGHT PATH ANGLE) HAS CONVERGED 

90 THETA* THETA+57 .2957 
GAMMA* GAMMA«57. 2957 
PSID* PSI*57.2958 
RETURN 
END 



SAMPLE CASE: VARIATION OF WIND-PROFILE MAGNITUDE 


aircraft aero. 

PROPFLLER 


POWER 


WINDS 

START POINT 

VARIABLE 

SET 

CODE 

W/S» lAA.O N/M2 

TS- 172.0 *^/S 


P/S* 

l.iO KW/M2 

WK« 0.00 

■ 

X 

o 

o 

o 

KM 

FIRST- 

0. 000 

N1 - 3 

N2- 10 

CL* .90 

S/A-p. 2.653 


RR» 

50.0 KM 

r^.U* 90.0 OEG 

HS- 

1 8.00 

KM 

FINAL* 

1.000 

N3- 50 
N4 ■ 1 

L/D* 36.6 


MAX 

P/W* 

6.62 KW/KN 



25.00 

KM 

STEP- 

.200 


L/D(ri)» .AIB 


MIN 

P / if^* 

.25 X MAX P/W 









WK 

X 

H 

R/C 

P/W-P 

P/W--S 

GAMMA 

THETA 

R 

VG 

VT 

VEC 

T 

AK 

ETA 

CP 

J 

P5I 


KM 

KM 

M/S 

W/N 

W/N 

DEG 

DEG 

KM 

M/S 

M/ S 

M/S 

SEC 





DEG 


0.00000 

0.00 

1 6.00 

.76 

2.73 

0.00 

.85 

155.6 

43.86 

51.3 

51 .3 

16.2 

0. 

.021 • 

716 

.021 

.937 

0.0 

0.00000 

. 51 

18 .01 

.78 

2.73 

0.00 

• 85 

155.8 

43.86 

51.3 

51.3 

16.2 

10. 

.021 . 

716 

.021 

.937 

0.0 

C. 00000 

27.52 

19.69 

6.18 

P .62 

1.95 

6.01 

123.4 

23.75 

59.0 

59.0 

16. 1 

510. 

.029 . 

889 

• 086 

1.078 

0.0 

0.00000 

60.98 

22.79 

4.16 

6.58 

0.00 

3.19 

A8. 3 

30.44 

74 .6 

74.6 

16.1 

1010. 

.045 . 

907 

.105 

1.363 

0.0 

0.00000 

8 3.61 

23 .41 

.67 

2.66 

0.00 

.54 

28.3 

49.43 

78 .7 

78.7 

16.2 

1312. 

.051 . 

902 

.047 

1.438 

0.0 

0.00000 

84.32 

23.41 

. 60 

2.60 

0.00 

.48 

27.9 

50 .05 

78.8 

78.8 

16.2 

1321. 

.051 . 

8 97 

• 046 

1.439 

0.0 

0.00000 

157. 9C 

21.68 

-1 .62 



-1.35 



6 6.6 

68.3 

16. 2 

2321 . 

-.0394 





0.00000 

222. 49 

20.16 

-1.44 



-1. 36 



60.9 

60.9 

16.2 

3321. 

-.0301 





0. 00000 

280.07 

16.79 

-i.3n 



-1.36 



54.6 

54.6 

16.2 

4321 . 

-.0239 





0.00000 

309.26 

16.09 

-1.24 



-i.37 



51.7 

51.7 

16. 2 

4871 . 

-.0211 





0.00000 

313. 19 

18 .00 

-1.23 



-1.37 



51.3 

51.3 

16.2 

4947. 

-.0212 






.20003 

0.00 

18.00 

.59 

2.60 

0.00 

.65 

155.8 

43 .86 

^0.2 

51.3 

16.2 

0. 

.022 . 681 

.020 


.937 

11.8 

.20000 

. 50 

18 .01 

. 59 

2.60 

0.00 

.65 

155.8 

43 .66 

50.2 

51.3 

16.2 

10. 

.022 .681 

.020 


.937 

11.8 

.2 0000 

26.93 

19.78 

6.19 

9.62 

1. 51 

6.03 

12^.7 

23.98 

57.9 

58.5 

16.1 

510. 

.029 .890 

.085 

1 

.068 

8.4 

.20000 

59.85 

22.70 

4. 37 

6.83 

0.00 

3.37 

49.8 

29.64 

73.6 

74.1 

16. 1 

1010 . 

.045 .904 

.107 

1 

.353 

6.4 

.20000 

83. 6A 

23.38 

.63 

2 .63 

0.00 

.51 

28.2 

49.44 

76.1 

78.5 

16.2 

1330. 

.050 .898 

. 0 46 

1 

.435 

6.3 

• 20C00 

84.34 

23. 39 

.57 

2.57 

0. 00 

.45 

27.8 

50 . 06 

78.1 

78 .6 

16.2 

1339. 

.050 . 893 

.045 

1 

.435 

6.3 

.20000 

157.30 

21 .66 

-1.6? 



-1.34 



6 8.2 

63.7 

16.2 

2339. 

0393 





. 20000 

221.25 

20.14 

-1. 44 



-1.36 



60.2 

60 • 8 

16.2 

3339. 

-.0300 





.20000 

278. OA 

id .77 

-1.30 



-i.37 



53.7 

54.5 

16.2 

4339. 

-.0250 





.20000 

'305. 7b 

13 .09 

-1.24 



-1.37 



50.6 

51.7 

16.2 

4871. 

-.0224 





. 200C0 

309 .69 

18.00 

-1.23 



-i.37 



50.2 

51 .3 

16.2 

4949. 

-.0225 







.^0000 

0.00 

13 .00 

.11 

2.21 

0.00 

. 12 

155.0 

43 .86 

. 4000 j 

.47 

16.00 

.11 

2.21 

0 . JO 

.12 

155. 

43. 6 6 

•40000 

24.^6 

19.37 

6.2 3 

6.62 

.01 

0.31 

123.9 ' 

24.~81 

•40000 

bo. 01 

22.36 

5.13 

7. 75 

0. 00 

4.07 

55.5 

27.06 

. 4COOO 

S3. 02 

23.26 

. 60 

2.59 

0.00 

.45 

26.4 

48.84 

•4U000 

84.39 

23.27 

.47 

2.47 

0.00 

.39 

27.7 

50.05 

•4COOO 

155.32 

21.56 

-1.61 



-1. 34 



.40000 

217. 19 

20.05 

-1.43 



-1.35 



•40000 

271. 40 

18.69 

-1 .30 



-1.36 



.40000 

293.92 

18.10 

-1. 24 



-1.38 



.40000 

297. 87 

18. OC 

-1 .24 



-1.38 




. 60000 

0.00 

16.00 

-1.29 

o.co 

1.55 

-i.39 

i55.6 

43.86 

.bOOOO 

. 40 

17.99 

-1.24 

0. 00 

1.55 

-1.39 

15 5.8 

43 .66 

. bOOOO 

20. OC 

18 .11 

2 .81 

4.60 

0.00 

3.11 

138.5 

27.27 

.60000 

44. SC 

20 .84 

6.08 

8.62 

2.23 

5. 52 

78.5 

21.20 

.60000 

77. 76 

22.69 

1 .00 

2.96 

0.00 

.77 

31.5 

43.45 

.60000 

34.66 

22.75 

.30 

2.27 

0.00 

. 26 

27.0 

50.06 

.60000 

150.21 

21.10 

-1.55 



-1.35 



.60000 

206.85 

19.64 

-1. 40 



-1.37 



.60000 

2 5 4. 42 

16 .31 

-1.27 



-1. 38 



• 60000 

262. 37 

18.07 

-1.25 



-1.36 



.oOOOO 

264. 57 

18.00 

-1.24 



-1.36 




. 80000 

0.00 

18. 

00 

-1.30 

0.00 

.82 

-i .^0 

155.8 

43 .86 

•80000 

.29 

17. 

99 

-1.25 

0.00 

.82 

-i. AO 

155.8 

43.86 

. 80000 

12. 50 

1 17 . 

37 

-1.20 

0.00 

.68 

-1. AO 

147.9 

32 .69 

WIND 

SPEED 

TOO 

LARGE 

AT 

16.9 

KN. 





UIND 

SPEED 

TOO 

LARGE 

AT 

16.9 

Kfi. 





WIND 

SPEED 

TOO 

LARGE 

AT 

16.9 

KM. 





WIND 

SPEED 

TDD 

LARGE 

AT 

16.9 

KM. 





WIND 

SPEED 

TOO 

LARG E 

AT 

16.9 

KM. 





W IND 

SPEED 

TOO 

LARGE 

AT 

16.9 

KM. 





WIND 

j P t E D 

TOO 

L aKG F 

AT 

16.8 

KM. 





WIND 

SPEED 

TOC 

LARGE 

AT 

16.8 

KM. 





WIND 

SPEED 

TOO 

LARGE 

AT 

16.8 

KM. 





WIND 

SJ^EED 

TOD 

LARGE 

AT 

16.9 

KM. 






WIND SPEED 

TOO 

LARGE AT 

16.0 

KM. 




W K 

XC 

H” 

TC 

E/W-S 

P/W-T 

XT 

TT 


KM 

KM 

HR 

KJ/N 

KJ/N 

KM 

HP 

0 . o 0 0 

84.32 

23 .41 

.367 

1.452 

9.803 

313.19 

1. 37 

.200 

84. 34 

2 3.38 

.37? 

1.377 

9.730 

309.69 

1.37 

• 4 C 0 

84. 39 

23.27 

• 3 8 9 

1.15 4 

9. 503 

297.87 

1.37 

• j ■»- J 

8m. 6 1 

22 .7 5 

.4^9 

1.1^7 

9.203 

264.57 

1.35 

.■200 

— i. 

— J 

T 

1 

- I 

-I 

17.38 

1. 35 

i .000 

_ 7 

-1 

- I 

-i 

-I 

0.00 

1. 35 


46.8 

51.3 

16.2 

0 . 

. 026 . 579 

.017 

.937 

24.2 

46.8 

51.3 

16.2 

10 . 

.026 .579 

.017 

.937 

24.2 

iz.y 

* 56 . 6 

16 . f ‘ 

sTo. 

.031 .891 

' . 079 

1.034 

18.5 

70.2 

72 . 1 

16.1 

1010 . 

.042 .892 

.115 

1.317 

13.1 

76.0 

77.8 

16.2 

1384 . 

.046 .891 

.045 

1.422 

12.5 

76.0 

77 . 9 

16.2 

1402 . 

.046 .882 

.043 

1.423 

12 . 5 

66.1 

68.1 

16.2 

2402 . 

-.0386 




58.1 

60 . 4 

16.2 

3402 . 

-.0300 




50.6 

54.2 

16.2 

4402 . 

-.0281 




47.3 

51 .7 

16.2 

4862 . 

-.0261 




46.8 

51.3 

16.2 

4946 . 

-.0262 





40.4 

51.3 

16.2 

0 . 

-.032 

0.000 

0.000 

0.000 

38.0 

40.4 

51.3 

16.2 

10 . 

-.032 

0.000 

0.000 

0.000 

38.0 

41.0 

51.6 

16.2 

510 . 

.032 

.063 

.035 

.942 

37.4 

56.7 

63.7 

16.1 

1010 . 

.034 

.891 

.100 

1 . 163 

22.7 

70.1 

74.3 

16 . 2 

1510 . 

.045 

.900 

.047 

1.357 

19.3 

70 . o 

74.7 

16.2 

1617 . 

.045 

.852 

.036 

1 . 365 

19.2 

60.9 

65 • 7 

16.2 

2617 . 

-.0356 




52.4 

58.4 

16.2 

3617 . 

-.0370 




42.9 

52 . 6 

16 .2 
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APPENDIX B 


GLIDE-TIME PARAMETER 


An expression for the time required to glide between two altitudes is given 
as equation 29 of reference 8. The development of that equation assumes, that 
the aerodynamic characteristics (Ci and Cn) remain constant and that accelera- 
tion effects (eq. (4)) are negligible. Tnat endurance equation for gliding 
flight is written as: 


L / L 
“ dVw/S (cos )-3/2 

- "1 ■ 

where h. and ho are the final and initial altitudes, respectively. 

Equation (Bl) can be simplified in two ways. First, since y is a small 
angle, the cosine term can be approximated as 1.0. Second, if the range of 
altitudes lies between about 16 and 26 kilomters, equation (11) can be used to 
approximate density variation by choosing a = 0.105 and b = 0.0013 
throughout that altitude range. 

Substituting equation (11) into equation (Bl) yields an integrable 
expression: 


■02 



(Bl) 


tg = e 

^ D V 2 W/S 


(a2/8b) 


*h2 


-(b/2)(h + (a/2b))2 

e 


^ hi 


dh 


L 


D 


1 (a2/8b) 




dz 



, (a2/8b) [~ 


(B2) 


(B3) 


where z = /b/^ (h + (a/2b)) 

Equation (B3) may be rearranged to produce an expression independent of vehicle 
aerodynamic characteristics. After substituting the values of a and b, the 
equation becomes: 

tgVgi = 27.873 (erf (zg) - erf (z^)) (B4) 
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where z * 1.0296 + 0.025495h 


where h Is expressed in kilometers, Vg in meters per second, and tq in 

hours. As in equation (Bl), h, is the final altitude because of the negative 
rate of climb. 

Glide time can be determined for a specific vehicle where L/D and Mq are 
given. For the class of vehicles considered in this study, the values of 
L/(DVe) lie approximately between 10 and 0.1. The largest value yields the 
longest glide time and is produced by low W/S and high L/D. 
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TABLE I.- DESCRIPTION OF BASELINE CONFIGURATION OF 


HIGH-ALTITUDE AIRPLANE PLATFORM 

Airplane aerodynamics: 

Aspect ratio, A 30 

Lift coefficient, C^ q ^9 

Lift-to-drag ratio 

Altitude function, L/D 36.6 + 0.418 h 

Folded propeller decrement, L/D 1.5 

Oswald efficiency factor, e 0.96 » 

Wing loading, W/S 144 W/vr 

Propel lers(s): 

Activity factor 80. 

Design lift coefficient 0.3 

Ratio of wing area to propeller-disk area, S/k 2.653 

Tip speed, V^ip i 172. m/s 

Motors (s) : 

Maximum specific power (available), 8.62 W/N 

Minimum specific power (required), (P/W)^^^,, 2.16 W/N 

Power transmission: 

p 

Power intensity at reference range, P/S 1.10 kW/m 

Reference range, R 50 km 

Range-power attenuation factor R/r 

Transmission initiation point 

Altitude, h 18 km 

Horizontal range, x 40 km 

Transmiss ion-termination^ si ant range 50 km 
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TABLE II.- INPUT PARAMETERS FOR PERFORMANCE PROGRAM OF APPENDIX A 


Array number 

Program name 

Description 

1 

WOS 

W/S 

2 

Cl 

Cl 

3 

BLOD 

L/D component independent of altitude 
of altitude 

4 

HLOD 

coefficient of altitude-dependent 
term in L/D equation, km"l 

5 

TS 

Vtip 

6 

SOAP 

S/Ap 

7 

POS 

P/S 

8 

RR 

R 

9 

POWL 

maximum P/W used by propulsion system 

10 

WK 

kw 

11 

AMU 

y 

12 

XS 

Xs 

13 

HS 

hs 

14 

POWR 

ratio of minimum P/W to maximum P/W 
for propulsion system 

15 

TDOT 

Td/T 

16 

HI 

hi 

17 

RKR 

kr 


N1 

code for flight mode calculation 
(1- single, 2- climb, 3- climb plus 
glide, or, 4- glide only) 


N2 

element in input array to be varied 


N3 

number of calculation points per 
listed line 


N4 

transmitter polarization code 
(1- linear; 2- circular) 


SF 

final value of variation set 


SI 

initial value of variation set 


SS 

step size of variation set 


38 










TABLE III.- OUTPUT PARAMETERS 


Output listing name* 

X 

H 

R/C 

P/W-P 

P/W-S 

GAMMA 

THETA 

R 

VG 

VT 

VE 

T 

AK 

ETA 

CP 

J 

PSI 

XC 

HC 

TC 

E/W-S 

E/W-T 

XT 

TT 


given in listing sequence 


PERFORMANCE PROGRAM OF APPENDIX A 










Ground track for perfect 
alignment of transmitting 
antenna with rectenna — . 



(c) Velocities in lateral plane. 
Figure 1. - Concluded 
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99% profile (reference 5) 
95% profile (reference 5) 
99% profile (reference 19) 
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10 m/s (19 kt) 
15 m/s (29 kt) 
20 m/s (39 kt) 


wind profile not exceeded 
99 % of the tim e (ref. S) - 


(a) True airspeed profile 


Figure 7. - Comparison of limiting wind profile and vehicle airspeed 
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Figure 8. - Effect of wing loading and lift-to-drag ratio on the 
performance of the baseline configuration. 
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17. - Climb performance for three values of k ; baseline 
configuration; x. = 40 km. 
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Figure 19. - Effect of variation in reference range and motor size. 
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Figure 20. - Effect of initiating climb at different altitudes and 
horizontal ranges. 
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(b) Mean and standard deviation profiles for April 


Figure 26. - Concluded. 
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(a) Variation of wind-profile magnitude 

Figure 28. - Effect of wind direction and profile amplitude on the 
performance of the baseline configuration. 
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